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ABSTRACT 

We investigate the structure and kinematics of the circumstellar disk around the Herbig Ae star 
MWC 758 using high resolution observations of the 12 CO (3-2) and dust continuum emission at the 
wavelengths of 0.87 and 3.3 mm. We find that the dust emission peaks at an orbital radius of about 
100 AU, while the CO intensity has a central peak coincident with the position of the star. The 
CO emission is in agreement with a disk in keplerian rotation around a 2.0 M star, confirming that 
MWC 758 is indeed an intermediate mass star. By comparing the observation with theoretical disk 
models, we derive that the disk surface density E(r) steeply increases from 40 to 100 AU, and decreases 
exponentially outward. Within 40 AU, the disk has to be optically thin in the continuum emission 
at millimeter wavelengths to explain the observed dust morphology, though our observations lack the 
angular resolution and sensitivity required to constrain the surface density on these spatial scales. 
The surface density distribution in MWC 758 disk is similar to that of "transition" disks, though 
no disk clearing has been previously inferred from the analysis of the spectral energy distribution 
(SED). Moreover, the asymmetries observed in the dust and CO emission suggest that the disk may 
be gravitationally perturbed by a low mass companion orbiting within a radius of 30 AU. Our results 
emphasize that SEDs alone do not provide a complete picture of disk structure and that high resolution 
millimeter- wave images are essential to reveal the structure of the cool disk mid plane. 



1. INTRODUCTION 

Millimeter-wave observations of circumstellar disks 
provide the best tool to measure the distribution and 
kinematics of the circumstellar material. These ob- 
servations have for years achieved 1-3" angular reso- 
lution to resolve the disk emission on scales down to 
an orbital radius of about ~70 AU at the distance of 
the n earby star forming reg ions ( e.g. [K oerner & Sa rgentl 
119951 : iGuilloteau fe Dutrevl 119981 : Ikitamura et al.ll2002D . 
More recent observations are reaching subarcsecond res- 
olution that are resolving disks on scales of the Kuiper 
Belt and yielding unprecedented informations on the 
innermost disk re gions (see, e.g, lAnd rews et al.l 
IBrown et all 1 2009: Hughes et al.~2009; Isel la et al.l 
120101 : iPietu et al.ll2006D . 

In many cases millimeter observations confirm the ra- 
dial distribution of dust and gas inferred from the anal- 
ysis of the disk spectral energy distribution (hereafter 
SED). The large majority of pre- main sequence stars 
exhibit a flux excess over the stellar photosphere from 
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near-infrared to centimeter wavelengths. This excess 
arises from a "classical", or primordial, disk that ex- 
tends from few stellar radii to few hundreds of astronom- 
ical units. However, ab out 10% of the observed systems 
(|Muzerolle et al.l l2010f ) exhibit a deficit of flux in the 
near- and mid-infrared compared to a "classical" disk, 
but have simi lar level of far-inf rared emission (i.e., "tran- 
sition" disk: IStrom et al.l lT989). The SED shape suggests 
that these stars lack the emission from warm dust close 
to the central star. The lack of warm dust emission may 
result from several different effects such as the presence of 
a stellar companion, the formation of a planetary system, 
disk photoevaporation, magneto -rotational instabilities , 
disk vi scosity, and grain growth ([Alexander & Arm itage 
20071: iChiang fc Murrav- Clavl 120071: llreland fc Kraus 



2008; llsella et al.l 120091 : iStrom et all 119891 : iTanaka et a! 
2005). Subsequent high-resolution millimeter- wave im- 
ages have confirmed that several "transition" disks are 
indeed dust depleted wi thin the radius predicted by 
the SED modeling E^iLeLalJl2QQ5 IBrown et al1l2009t 
iHughes et alll2009h . 

The increasing sensitivity and resolution of 
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(sub) millimeter- wave interferometers is now pro- 
viding unexpected insights about the structure of 
"classical" disks. Sub-arcsecond observations of AB Aur 
and RY Tau systems show indeed that these two 
circumstellar disks may be dust deple ted within orbital 
radii of 70 AU an d 13 AU respectively (jPietu et al.ll2Q05b 
Ilsella et al.l 12010). although no disk cl earing has been 
inferr ed from the analy sis of the SED (jRobitaille et al.l 
120071 : iDullemond et al.l 1200 if ). It is therefore possible 
that other "classical" disks may be characterized by 
a disk surface density that drastically deviates from 
that inferred from the SED modeling. In particular, 
since the dust opacity in the infrared is much larger 
than at millimeter- wavelengths, it is reasonable to 
assume that a large fraction of "classical" disks might by 
characterized by partially cleared inner cavities which 
are not detectable through the analysis of the infrared 
SED. 

In this paper we discuss the Herbig Ae star MWC 758 
(HD 36112), a 3 Myr old pre-main sequence A5 star 
located near the edge of the Taurus star forming re- 
gion at a distance of 2001^ P c (|Perrvman et al.l 1997; 
Ivan den Ancker et al.l Il998f ). Previous investigations 
based on the SED modeling and on spatially resolved 
interferometric observations suggested that MWC 758 
is surrounded by a classical, fully fl ared rotating disk 
that extends from few stellar radii (Isel la et al.l 120081 : 
Eisn er et al.l 120041) up to an outer radius of 385 ± 26 
Alfl (Chapillon et al. 2008). At 1.3 mm, dust continuum 
emission observed at 2" angular resolution is centrally 
peak and more comp act than the 12 CO (2-1) emission 
(jChapillon et al.ll2008f ). This difference in radial exten- 
sion between dust and gas emission is sim ilar to that 
found in other circum stellar disks (see, e.g.. iPietu et al.l 
[20051 : Ilsella et al.l f2007h and can be explained with an 
exponential decrease of the disk surface density in the 
outer most disk regions (Hug hes~et al.l 120081 : Ilsella et all 
l2010h . 

Here we present new images of the millimeter-wave 
dust continuum and 12 CO emission toward MWC 758 
characterized by an angular resolution of 0.7", which 
is three times better than previous observations. Our 
observations resolve the inner disk into two peaks in 
the continuum emission, strongly suggesting that the 
disk surface density decreases within an orbital radius 
of about 100 AU. Similarly to the cases of AB Aur and 
RY Tau, the MWC 758 inner disk may be dust de- 
pleted with respect to a power law surface density profile, 
though no disk clearing has been inferred from the anal- 
ysis of the SED. 

The observations, presented in Section [2] and [3j are 
compared with theoretical disk models to constrain the 
surface density profile of MWC 758 disk. The results are 
presented in Section HI The discussion and conclusions 
follow in Section [5] 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. SMA observations 




1 1 1 1 1 1 1 1 1 1 

Obs (3.3 mm) 




I 1 1 1 1 I 1 1 1 1 
Mod (3.3 mm) 





1 1 1 

- Obs 


1 1 1 1 I 1 1 1 1 1 1 1 1 
- Mod (3.3 mm) 




(g) 











1 -1 

6R.A. (arcsec) 

Figure 1. Observations, best fitting models and residual maps 
of the dust continuum emission observed toward MWC 758 by 
SMA at 0.87 mm (left column) and by CARMA at 3.3 mm (right 
colums). The la noise levels are 2.1 and 0.25 mjy/beam in SMA 
and CARMA observations respectively. Contours are plotted every 
3cr. The symbol * marks the position of MWC 758 after correcting 
for the proper motion as described in Section [3] The dashed line 
shows the orientation of the disk major axis derived from the CO 
observations as discussed in Section [3] 

The MWC 758 system was observed at 345 GHz with 
the Submillimeter Arra}0 (SMA) on Mauna Kea, Hawaii, 
on 2008 Jan 28. The array was in the extended configura- 
tion, and the seven available antennas provided baseline 
lengths from 31 to 225 meters. The atmospheric trans- 
parency was excellent through the track, with an opacity 
of 0.04 at 225 GHz measured at the nearby Caltech Sub- 
millimeter Observatory. The correlator provided 2 GHz 
of bandwidth in two sidebands separated by 10 GHz, 
with the upper sideband centered near the CO J=3-2 line 
at 345.796 GHz. This line was observed at a spectral res- 
olution 0.18 km s _1 in a single correlator "chunk" span- 
ning 104 MHz. The primary flux calibrator was Titan, 
and the passband calibrators were the quasars 3C273 and 
3C111. Complex gain calibration was carried out with 
observations of the quasar J0530+135 (1.75 Jy) inter- 
leaved with MWC 758. Data was obtained for MWC 758 



1 Note that Chapillon et al. and Eisner et al. assume a stellar 
distance of 140 and 150 pc respectively. The radii quoted from 
these two papers have been therefore rescaled for the distance of 
200 pc 



2 The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics, and is funded by the Smith- 
sonian Institution and the Academia Sinica. 
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over the hour angle range —3.2 to +4.7. All of the data 
editing and calibration were performed using standard 
routines in the MIR software package for IDL. The re- 
sulting map of the dust continuum emission is shown in 
Figure [TJ Applying natural weighting to the data we 
achieve a noise level of 2.1 mJy/beam and synthesized 
beam FWHM of 0.76"x0.56". Figure [2] shows the CO 
emission in the velocity range between 2.94 and 8.93 km 
s _1 , as obtained by averaging the observations every two 
channels. The resulting noise level is 83 mJy/beam. 

2.2. CARMA observations 

The dust thermal emission toward MWC 758 was ob- 
served on January 2009 using the Combined Array for 
Research in Millimeter-wave Astronomy (CARMA) lo- 
cated at the altitude of 7200 feet in the Inyo mountains 
of eastern California. The array was in the extended 
B configuration, which provides baselines between 80.9 
and 946.3 m. The observations were obtained at the 
local oscillator frequency of 97.5 GHz (A ~ 3.3 mm). 
The CARMA correlator was configured with six bands, 
each of which was configured to 468 MHz band-width 
to provide maximum continuum sensitivity. No molecu- 
lar emission line was observed with CARMA during this 
track. The passband shape was calibrated by observing 
3C454.3 and the flux calibration was set by observing 
Uranus. The quasar J0510+180 was observed every 15 
minutes to correct for atmospheric and instrumental ef- 
fects. The data reduction and the image reconstruction 
were obtained using the MIRIAD software package. Fig- 
ured] shows the resulting map, which is characterized by 
a rms noise of 0.25 mJy/beam and a synthesized beam 
FWHM of 0.80" x 0.62" when natural weighting is applied 
to the data. 

3. RESULTS 

Figures [TJ [2j and [3] show the dust and CO emission 
towards MWC 758 as observed by pointing the telescope 
at the ICRS coordinates of the star at the epoch J2000, 
RA=05h 30m 27s.530 (±12.89 mas) and Dec=+25° 19' 
57". 08 (±5.57 mas). The position of the star at the date 
of the observation is shown by the symbol * and is calcu- 
lated by applying a proper motion of c)(R.A)=5.23±1.40 
mas/yr and £(Dec)=-25.95±0.66 mas/yr to the J2000 
coordinates (bot h the stellar coordinat es and the proper 
motion are from iPerrvman et al.lll997f ). The uncertain- 
ties on the stellar position resulting both from the Hyp- 
parcos astrometry and proper motion uncertainties are 
smaller than the size of symbol itself. 

The dust emission at both 0.87 and 3.3 mm is charac- 
terized by two peaks which are symmetric with respect of 
the position of the star and are separated by about 0.7", 
or 140 AU at the distance of 200 pc. At the wavelength 
of 0.87 mm, the northern and southern peaks have fluxes 
of 89 and 77 mJy/beam respectively. A similar flux ratio 

3 Support for CARMA construction was derived from the Gor- 
don and Betty Moore Foundation, the Kenneth T. and Eileen L. 
Norris Foundation, the James S. McDonnell Foundation, the As- 
sociates of the California Institute of Technology, the University 
of Chicago, the states of California, Illinois, and Maryland, and 
the National Science Foundation. Ongoing CARMA development 
and operations are supported by the National Science Foundation 
under a cooperative agreement, and by the CARMA partner uni- 
versities. 



between the two peaks is measured at 3.3 mm, suggesting 
that the observed morphology corresponds to an asym- 
metric dust emission between the two sides of the disk. 
By integrating over a circle with a radius of 1.5" centered 
at the position of the star, we measure fluxes of 217±40 
mJy and 4.9±1.7 mJy at 0.87 and 3.3 mm respectively, 
where the errors account for both the thermal noise and 
absolute calibration uncertainties. The integrated fluxes 
are in agreement with the values meas ured at 2" angu- 
lar resolution by iChapillon et al.l (|2008f ) suggesting that 
the dust emission is all coming from the spatial scales 
probed by our observations. Assuming that the disk mil- 
limeter emission follow a power law F(v) oc v a , we derive 
a spectral index ai_3 mm =3.0. 

Figure [2] shows the 12 CO J=3-2 channel maps in the 
velocity range between 2.94 and 8.93 km/sec. The CO 
emission is spatially resolved and shows a velocity pat- 
tern typical of an inclined rotating disk. The first mo- 
ment of the CO emission, i.e., the velocity gradients 
across the disk, indicates that the western and eastern 
sides of the disk are respectively red- and blue-shifted 
with respect to the relative velocity of the star (see the 
left panel of Figured]). The velocity gradient also sug- 
gests a position angle of the disk major axis of about 
65° as measured from north to east. The central panel 
of Figure [3] shows the spectrally integrated CO emission. 
In contrast to what is observed in the dust continuum, 
the CO emission has a single peak coincident with the 
position of the star. The CO emission is observed up 
to a distance of 1.3" in the red-shifted side of the disk, 
and only to about 1" in the other side. The line in- 
tegrated flux is 12.9 Jy km/sec. This value is about 
14% sm aller than that m easured in single dish obser- 
vations (|Dent et al.l [2005). suggesting that part of the 
CO emission is resolved out in the inter ferometric ob- 
servations. As a consequence, the observed asymmetry 
in the CO emission might be due to the lack of short 
baselines in the SMA observations instead of to a real 
asymmetry in the gas temperature or distribution. Fi- 
nally, the spatially integrated CO spectrum (right panel 
of Figure [3]) has a double peak profile centered at the ve- 
locity of vi sr = 6.0 km/sec. The line emission is observed 
in the velocity range vi sr ±4 km/sec with a brighter flux 
in the red-shifted half of the disk. 

To date, the different morphologies observed in the 
dust and gas emission toward MWC 758 represent a 
unique case among the sample of circumstellar disks spa- 
tially resolved with millimeter- wave interferometers. Our 
observations suggest that the optical depth of the dust 
continuum emission decreases within a radius of 70-100 
AU from the star. Despite that, the disk remains gas- 
rich down to at least a distance of 30 AU to explain the 
centrally peaked CO emission. In addition, the SED of 
MWC 758 shows a flux excess above the stellar photo- 
sphere between near- and far-infrared wavelengths typi- 
cal of optically thick flaring disk extending to few stellar 
radii (see Figure [5]). The next section is devoted to recon- 
ciling all these observations in the framework of existing 
disk models, in order to constrain the radial distribution 
of the circumstellar material around MWC 758. 

4. DATA ANALYSIS 
4.1. Disk models 
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Figure 2. Channel maps for the 12 CO J=3-2 line emission observed toward MWC 758 in the velocity range between 2.94 and 8.93 km/sec. 
The contours are spaced by 3a = 0.25 Jy/beam. The figure also shows the residuals obtained by subtracting the best fitting model CO 
emission from the observations. The symbol * marks the position of the central star after correcting for the proper motion. 
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Figure 3. Left panel: 12 CO J=3-2 velocity gradient across MWC 758 circumstellar disk. Red and blue colors indicate the red- and 
blue-shifted disk velocities respectively. Velocity contours start at 4.7 km sec -1 and are spaced by 0.3 km sec -1 . The dashed line indicates 
the position angle of the apparent disk major axis, which is rotated by 65° to east from north. Central panel 12 CO J=3-2 integrated 
intensity map. Right panel: 12 CO J=3-2 spectrum obtained integrating over the region shown in the intensity map. 



To analyze the data we a dopt the disk model discussed 
in llsella et"aH (|2007L 12009). In brief, we compute the ra- 
dial dust temper ature profile following the "two-layer" 
approximation of lChiang fc Goldreichl (|1997[ ) , where the 
disk is characterized by a warm surface layer and a cooler 
interior. Both temperatures are calculated as a func- 
tion of the or bital radius by iterating on the vertical disk 
structure (see lDullemond et al.ll2001f ). The disk is in hy- 
drostatic equilibrium between the gas pressure and the 
stellar gravity, which leads to a flared geometry with the 
opening angle increasing with the distance from the cen- 
tral star. 

The dust opacity is c alculated by assumi ng an interstel- 
lar grain composition (Pollack et al.lH994j ) and a particle 
size distribution between a minimum and a maximum 
value dmin and a max according to n(a) oc a~ q . We fix 



dmin — 0-05 /im and vary a max and q to reproduce the 
spectral index a of the millimeter disk emission. For sake 
of simplicity we assume that the d ust opacity is constan t 
throughout the disk (although see iBirnstiel et al.1 l2010). 

The radial distribution of the circumstellar material 
follows the similarity soluti on for the disk surface density 
of a viscous keplerian disk ([Lvnden-Bell fc PringL 
expressed by 



e density 

na [T97i 



E(r,t) 



xexp 



1 



n 



(2-7) 



1 



2(2- 7 ) 

(1) 

where the characteristic radius r*t, 7, and the surface den- 
sity normalization E t are free parameters of the model, 
as well as the disk inclination and position angle. 
From the derived dust density, temperature, and opac- 
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Table 1 

Best fitting model parameters for the dust and CO millimeter 
emission observed towards MWC 758. The quoted uncertainties 
correspond to the la uncertainty level 

Disk orientation: 

Inclination = 21° ±2°, P. A. = 65° ±7° 

Surface density H(R): 

7 = -4.1±0.3, R t = 69±3 AU, gt =1.5±0.1 g cm" 2 
CO temperature T(CO): 

T C q(50 AU) = 53±4 K, C = -0.1±0.2 

M* = 2.0±0.2 M 

Grain size distribution n(a): 

amax = 1 cm, q = 3.5 

ity we calculate the disk SED and synthetic disk images 
in the dust continuum at 0.87 mm and 3.3 mm using the 
radiat ive transfer solution discussed in iDullemond et al.l 
(2001). The synthetic disk images are then Fourier trans- 
formed and sampled at the appropriate positions on the 
(u,v) plane corresponding to our CARMA and SMA ob- 
servations. 

The CO emission is calculated assuming that rota- 
tional levels are thermalized. The CO temperature is 
parameterized as a function of the radius as Tco(r) = 

^co( r o)( r o/ r )^ where Tco(^o) an d C are f ree parame- 
ters and ro = 50 AU. The CO density is then calculated 
by assuming that the dust and the gas are well mixed 
with a gas-to-dust ratio of 100. The adopted CO abun- 
dance is 12 CO/H 2 = 7.0 x 10- 5 . Finally, the CO channel 
maps are calculated in the velocity range probed by our 
observations by assuming that the emission originates 
from an inclined keplerian disk. Both the disk inclina- 
tion and the mass of the central st ar are free parame ters 
of the model (see the Appendix of llsella et aLll2007l for 
more details on the CO emission model). As for the case 
of the dust emission, the synthetic CO channel maps are 
then Fourier transformed and compared to the measured 
complex visibilties. 

The best fitting models for both the dust and gas emis- 
sion, and the uncertainties on the model free parameters, 
are then obtained by y 2 minimization carried out us- 
ing the Markov C hain Monte Carlo method discussed in 
llsella et al.l (|2009D . 

4.2. Results of the model fitting 

A satisfactory fit of the spatially resolved dust and CO 
emission, as well as of the SED between far-infrared and 
(sub) millimeter wavelengths is obtained with the disk 
model presented in Table [TJ The disk surface density 
is mainly constrained by the optically thin dust emission 
at 0.87 and 3.3 mm and is characterized by 7 = —4.1 
and r t = 69 AU. As shown in Figure gj E(r) has a "bell" 
shape that reaches the maximum value of 3.3 g/cm 2 at 
the radius of 100 AU. For r < r t , the surface density is 
proportional to r 4 , suggesting that a steep radial density 
profile is required to explain the two peaks observed in 
the dust emission. For r > 100 AU, E(r) fades expo- 
nentially as prescribed by Equation [TJ Given the angu- 
lar resolution and sensitivity of our observations, we can 
constrain the surface density only between 40 and 130 
AU. In this range, the disk is optically thin to the dust 
continuum emission at both 0.87 and 3.3 mm. Within 40 
AU, the radial profile of S(r) might significantly deviate 
from that shown in Figure HJ However, the disk must 
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Figure 4. The solid red line shows the best fitting model solution 
for the disk surface density S(r) obtained by analyzing the spatially 
resolved observations of the dust and CO emission at both 0.87 and 
3.3 mm. The dashed lines indicate the 3cr uncertainty range. Given 
the angular resolution and sensitivity of our observations, S(r) is 
constrained only between 40 and 130 AU. Within 40 AU, S(r) can 
significantly deviate from the profile shown in the figure, though 
the disk must remain optically thin at both 0.87 and 3.3 mm to 
explain the two peaks observed in the dust emission. 

remain optically thin to the dust continuum millimeter- 
wave emission down to at least few AU in order to re- 
produce the two peaks shown in Figure [TJ 

Radially integrating the surface density we obtain a 
disk mass of 0.01 M . This result is affected by a large 
error due to the uncertainties on the dust opacity and 
on the gas to dust ratio. In particular, the analysis of 
the 12 CO (2-1) and (1-0) line emission suggests that the 
gas to du st ratio in MWC 758 disk might me smaller 
than 100 (jChapillon et al.l l2QQ8h . leading therefore to a 
smaller disk mass. 

Despite the steep surface density profile, the 12 CO (3- 
2) emission remains optically thick down to a radius 
of about 20 AU. As a consequence, the CO emission 
does not constrain the disk surface density but it pro- 
vides a measure of the CO temperature in the form 
Tco(r) ~ 53K x (SOAU/rT - 1 . Note that the 12 CO 
(3-2) emission remains optically thick even if the CO 
is largely depleted with respect to the standard abun- 
dance. In particular, a ssuming a CO abundance of 10 -6 
(|Chapillon et alJ f2Q08) we calculate a mean optical depth 
of about 10 3 . The almost radially constant profile of the 
CO temperature is in good agreement with the value 
( = 0.05 ± 0.2 deriv ed from the ana l ysis o f the 12 CO 
(2-1) emission line by Chapill on et al.l (j2008[ ). In the as- 
sumption that the disk is in keplerian rotation, we derive 
a disk inclination of 22° and a dynamical mass for the 
central star of 2.0 M , similar to the value of 2.3 M 
derived from the H- R diag ram using stellar evolutionary 
tracks (jManoi et al.1120061 ). 

The residuals obtained by subtracting the best fit- 
ting model to the observations reveal asymmetries in the 
MWC 758 disk. In the case of the dust emission (see 
Figure [1]), our disk model fails to reproduce the flux mea- 
sured in the northern side of the disk at both 0.87 and 
3.3 mm. From the intensity of the residuals we calcu- 
late that the dust density, opacity, and/or temperature 
in the northern side of the disk might be 20% higher than 
the southern side. Figure [2] shows the residuals for the 
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CO channel maps in the velocity range between 2.94 and 
6.11 km s _1 . Most of the significative structures appear 
in low- velocity channels of the red-shifted side of the disk 
between 6.11 and 7.11 km/sec. Since the CO emission is 
almost completely optically thick, this suggest that the 
CO temperature might be higher in this side, or that the 
CO extends to slightly larger distances from the central 
star. Note, however, that the CO asymmetry might also 
be due to the lack of short spatial frequency in SMA 
observations as discussed in Section [3l 

Figure [5] shows the comparison between the observed 
SED and the SED of the best fitting model to the mil- 
limeter data. At millimeter wavelengths, a good agree- 
ment with the observations is obtained with a grain size 
distribution slope q = 3.5 and a maximum grain size 
Umax = 1 cm. Fitting the dust opacity fc(A) between 
1.3 and 2.8 mm with a power law fc(A) ex A _/3 , we ob- 
tain (3 = 1.35, which is a factor of two higher than the 
mean value measured in T Tauri stars but in agreement 
with those found in other Herbig Ae stars (see, e.g., 
iRicci et al]l201ttlNatta et al.ll2QQ4f i. 

The good agreement between the model and the obser- 
vation extends down to a wavelengths of about 20 fim. 
This suggests that the disk has a flared structure in the 
region that emits most of the mid- and far-infrared emis- 
sion, namely between 20 and 110 AU. At shorter wave- 
lengths, the emission of the best fitting model is much 
lower than what observed toward MWC 758. This is due 
to the fact that for £(r) ex r 4 , the dust density within 
few AUs from the central star is lower than 10 -6 g cm -2 
and its near-infrared emission is negligible with respect 
to the stellar photospheric emission. The observed near- 
infrared excess therefore suggests the presence of mate- 
rial with temperature between 1000 and 2000 K in excess 
to what predicted by our model. 

To reconcile the near-infrared SED with our 
millimeter-wave observations we propose two different 
scenarios. In the first case, the near-infrared emission 
is coming from an optically thick disk extending be- 
tween 0.05 and 5 AU. This model comprises a gaseous 
disk between 0.05 and 0.1 AU, a "puffed- up" inner rim 
of the dusty disk b etween 0.4 and 0.5 AU (see, e.g, 
llsella fc Nattall2005h . and a flaring disk between 0.5 and 
5 AU. As discussed in llsella et al.l (|2008f ). this model ex- 
plains the spectro-interferometric observations of the H- 
(1.6 /iin) and K- band (2.2 /im) disk emission (see also 
iBenistv et all 1201(1 ). The resulting SED, which is indi- 
cated by the short-dashed green line in Figure [5j provides 
a good description of the near-infrared emission and of 
the 10 /am silicate feature. Alternatively, the infrared 
SED can be explained with the presence of an optically 
thin layer of dust surrou nding the star (jVinkovic et al.l 
2006; Mulders et al. 2010). In this case, we calculate that 
a ring of dust extending between 0.15 and 0.8 AU, and 
characterized by an optical depth for the stellar radiation 
along the direction perpendicular to the disk mid plane 
of 0.4 is required to reproduce the disk emission (see the 
long-dashed curve in Figure [5}. For the assumed dust 
opacity, this model requires a minimum surface density 
of 0.01 g cm 2 , or a total mass of 2x 10~ 9 M . Note that 
neither of these models will produce a significant dust 
and CO emission at millimeter wavelengths. 

To summarize, assuming the similarity solution for a 
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Figure 5. Spectral energy distribution of MWC 758. Pho- 
tometric clata__fR^_the literature are shown with empty 
squares (Eisner et al. 2004; Malfait et al. 1998; Beskrovnava et al. 
1999; Cutri et al. 2003; Beichman et al. 1988; Elia et al. 2005; 
Chapillon et al. 2008). The thick black line show the observed IRS 
spectrum from the Spitzer archive. The filled squares show our new 
measurements at 0.8 and 3.3 mm. Optic al and near-infrare d data 
were de-reddened adopting Ay = 0.2 ( van den Ancker et al. 1998). 
The dotted black line shows the stellar photosphere modeled using 
a black body emission with a temperature T e f f = 8200 K, L* = 22 
L0, and a distance of 200 pc ( van den Ancker et al. 1998). The 
solid red line shows the SED for the best fitting model discussed 
in Section [3] The long-dashed blue and short-dashed green curves 
show the SED o f an optically thick and thin inner disk model, as 
discussed in Sec. 14.21 

keplerian viscous disk, we find strong evidence that the 
surface density of the MWC 758 disk increases as r 4 
between 40 and 100 AU from the central star, and de- 
creases exponentially outwards. However, this solution 
for the surface density is not unique and observations 
with higher angular resolution are required to discrim- 
inate between different models. We note, for example, 
that a discontinuity in the surface density at about 80 
AU from the central star might provide an equally good 
fit with a different step gradient. Within 40 AU, S(r) 
is constrained by the analysis of the SED. We calculate 
that a minimum surface density value of 10 -2 g cm -2 
between 0.15 and 0.8 AU is required to reproduce the 
observed near- and mid-infrared emission. Probing these 
spatial scales requires high sensitivity 0.1" resolution ob- 
servations that will become possible with the Atacama 
Large Millimeter Array (ALMA). 

5. DISCUSSION AND CONCLUSIONS 

The surface density profile in MWC 758 is similar to 
that measured in the transitional disk around the T Tauri 
star LkCa 15. In this case £(r) peaks a t about 50 AU 
from the central star (jPietu et al.l 120061 : lEspaillat et al.l 
2007) and a minimum mass of 1 x 10 11 M within 5 AU 
from the central star is required to produce t he observed 
near-infrared emission (Muld ers et "all l2010f ). Similarly 
to what suggested for LkCal5, the shape of the surface 
density in MWC 758 may be due to the disk clearing 
operated by a low mass stellar companion, or by one 
or more giant planets. To this regard, we note that re- 
cent spectroastrometric observations of MWC 758 ex- 
clude the presence of a stellar companion more massive 
than about 0.9 M orbiting between 0.15" and 0.3", or 
30 and 60 AU for the ass umed stellar distance of 200 pc 
(Wheelwright et aT1l2010L Wheelwright H. private com- 
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munication). The nature of a possible companion of 
MWC 758 remains therefore uncertain. 

The asymmetries observed in the dust and CO emission 
support the hypothesis that the disk may be gravitation- 
ally perturbed by a massive body. The asymmetric CO 
emission can be due to a warped optically thick inner 
disk, i.e., more inclined than the outer disk, so that the 
CO in north-west side of the disk is partially shielded 
by the inner disk from the stellar radiation. This would 
result in a lower gas temperature and consequently lower 
emis sion with respect to the other side of the disk (see, 
e.g.. iPanic et al.lfeolOf ). Indeed, we note that the incli- 
nation of the inner disk derived from near-in frared inter- 
ferometric observations is of a bout 40° ±10° (jlsella et al.l 
2006, 120081 : lEisner et all 120041 ). which is larger than the 
21°±2° derived in this work. 

Other mechanisms, which are in principle able to evac- 
uate the inner disk, do not seem to be consistent with 
the observations. Photoevaporation and magneto rota- 
tional instability models predict completely cleared inner 
disks and surface density profiles at the inner disk ra- 
dius much steeper than r 4 (I Alexander fc Armit age 2007; 
iChiang fc Murrav-Clayl l2007f ). In addition, the surface 
density within few AU is much larger than that predicted 
by the viscous solution for the surface density that best 
fits the millimeter data. This fact excludes that S(r) 
is shaped only by the rad ial profile of the disk viscos- 
ity (see the discussion in llsella et aD l2QQ9h . Low dust 
opacity caused by grain growth has also been suggested 
to explain inner holes observed in the continuum emis- 
sion at millimeter wavelengths (Tana ka"et al.ll2005[ ). In 
the case of MWC 758, the dust opacity should decrease 
by at least two orders of magnitude to explain the ob- 
servations. This would require a grain size distribution 
flatter than what typically assumed (i.e., n(a) ex a~ q 
with q < 3) and a steep radial decrease of the maximum 
grain (i.e, form 10 cm at 10 AU to 10 /im at 100 AU). 

In conclusion, our observations show that the disk 
around MWC 758 is characterized by a steeply increasing 
surface density between 40 and 100 AU from the central 
star. Given the infrared SED and the asymmetries ob- 
served in the gas and dust emission, we conclude that 
the surface density might be shaped by the presence of 
stellar, or planetary size companion orbiting within 30 
AU. The fact that MWC 758 was previously labeled as 
"classical" disk emphasizes that SEDs alone do not pro- 
vide a complete picture of disk structure and that high 
resolution millimeter- wave images are essential to reveal 
the structure of the cool disk mid plane. 
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